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Abstract. An analytical result has been obtained for the value of the dynamical 
charge necessary for calculations of oxygen vacancy kinetics in dielectric perovskite- 
type crystals. It is shown by using the Berry phase analysis that this charge equals 
the nominal charge of the vacancy: for example, for the double charged state, it is 2; 
a neutral vacancy has the zero charge. 
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Oxygen vacancies in oxides of the perovskite family can be easily created or 
removed by thermal treatment in reduction or oxidizing atmosphere respectively 
The vacancies are always present in single crystals and even more so in ceramics and thin 
films. This nonstochiometry effects some properties of perovskite-type oxides: dielectric 
permittivity can have a noticeable 0, ^, ^ and sometimes even large ||^ contribution 
from electrons connected with the oxygen vacancies; Second Harmonic Generation and 
luminescence experiments revealed a large effect of oxygen vacancies on the intensities 
0] ; Electroconductivity of the samples is strongly dependent on the oxygen vacancy 
concentration an, hence, on thermal treatment [§. 

Sometimes the presence of oxygen vacancies helps to improve properties of 
materials, for instance, it is possible to obtain extremely high dielectric permittivity in a 
wide temperature interval due to Maxwell- Wagner relaxation It happens that these 
unavoidable point defects spoil some properties, for example, in ferroelectric memory 
applications the presence of oxygen vacancies is thought to be crucial and results in 
the phenomenon called "fatigue" |P, |T0|: ferroelectric thin films degrade after a large 
number of polarization switching due to the movement of the oxygen vacancies leading 
to the creation of large frozen internal fields preventing the switching process. 

In order to study the last effect one can use kinetic equations describing a charge 
transport in a bias field [|11|]. Charged oxygen vacancies are thought to be moved in this 
field 1^, |T0|. For quantitative estimations one needs to know the oxygen vacancy charge 
Z and the field Ea acting on the vacancy. In review |T^ this problem was discussed 
in details: there had been several trials to describe these quantities among which there 
are two very different proposals, one of them suggests using Lorentz's expression for the 
local field, Ea = {e + 2)E/3, and the other assumes that each of the vacancies can be 
surrounded by a sphere and the field inside this sphere is described by an Onsager's 
expression: Ea = 36/{2e+ 1)E where e is dielectric permittivity and E being the 
bias macroscopic field. Scott and Dawber [|l^ criticized the first approach as giving an 
unrealistic diverging value of the local field at large e inherent to ferroelectrics, and they 
suggested using the Onsager expression which much better suits experimental data. 

In this letter I will show that the problem of choosing the local field and charge of 
the vacancy can be solved analytically with the help of a Berry phase analysis [|l^]. This 
approach rigorously considers not only the local field effect but also the covalent effect 
connected with the covalent binding of the ions. The resultant polarization is computed 
quantum mechanically by taking integrals over the Brillouin zone. 

Instead of dealing with local fields it is easier to consider average, macroscopic, 
field. In this case the charge of the vacancy must be replaced by the dynamical charge 

dPa 

where Pa is a-th component of polarization, and r^^ is the displacement of z-th ion in 
the /3-th direction. As the average field is uniform, only the dynamical charge should 
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be averaged over a path of the oxygen vacancy 

Z=]J Z*y^,,{T)dr (2) 

where / is the length of this path; for the sake of simphcity, we directed the field along 
the z axis. By substituting definition ([1|) to the integral in (0) one has 

1 f dP, , AP, 
Z=j T^dr = 3 
I J orvoz Iz 

where AP is the finite difference in polarization after performing the run over the path, 
and Iz is the z component of the displacement. 

The difference in polarization AP is the sum of polarization stemming from the 
ionic charges transferred, APj = Zil^ where Zi is the nominal (ionic) charge of the i-ih 
particle, and the electronic polarization APg: 

AP, = APi + APe (4) 

The computation of the electronic (covalent) contribution can be carried out with 
the help of the Berry phase analysis : 

APe ~ Ayp (5) 

where if is the average over the electronic bands Berry phase. The Berry phase strongly 
changes inside the unit crystal cell and it is a periodic function with respect to the 
displacement with the period coinciding with the lattice parameter. 

Consider the path shown in Fig. 1. The initial and final states in this path have 
the same point symmetry and, consequently, the electronic contribution vanishes in this 
path: APe = 0. In other words, at the displacements by a lattice constant or by a few 
lattice constants, the electronic contribution to the polarization vanishes and only the 
ionic contribution remains. This rigorously proves that the average dynamical charge 
of the oxygen vacancy equals its nominal charge; for example, for the double charged 
vacancy Vj, this charge is 2, for the single charged vacancy it equals 1, and, finally, for 
the neutral vacancy it is zero. 

It follows from the analyses performed that the Lorentz expression for the local field 
provides wrong result for the kinetics of oxygen vacancies in oxides. This expression 
gives a very large value which, in fact, reflects the fact of strong local field in the 
cenrosymmetric position of the simple cubic lattice made of polarizable ions. If one 
averages this field over the unit cell, as it was evenly suggested in Ref. then the 



difference between the local field and the average one obviously vanishes at all. Notice 
that we have shown that it is enough even to average the result over any of the paths 
coinciding two points in the lattice with equal point symmetry. Hence the field averaged 
over such a path equals the field averaged over the unit cell volume. Close result can 
be obtained if one uses the Onsager expression at high e: this result is only 3/2 times 
larger with respect to the right expression which uses the average field instead of the 
local one. 
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Caption 

Fig. 1. A path of the oxygen vacancy between equal points in the perovskite-type 
lattice. Large circles correspond to B ions in the ABO3 perovskite structure, small ones 
to oxygens, and the square is the oxygen vacancy. 
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